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Abstract Group living has increasingly been emphasized 
due to its potential consequence on inbreeding, nevertheless, 
the relationship between group living and inbreeding risk 
is largely unknown. The endangered snake species, Shedao 
pit viper (Gloydius shedaoensis) inhabiting mostly on Shedao 
Island (meaning “snake island” in Chinese), provides an 
excellent model system for research on group living and 
inbreeding risk. Small island area, high population density 
and extreme seasonal foraging opportunity prompt many 
individuals to live on the same tree, which likely increase 
the potential mating among relatives. To confirm this 
probability, we used genotypes from 13 microsatellites 
DNA markers and examined the genetic relationships 
between pair of individuals lived on the same tree. The 
results showed that snakes on the same tree did not show 
closer relationships compared with individuals from 
different trees. The group constructions based on scenarios 
of parental-offspring, full-sibling, half-sibling as well as 
cousin were not consistent with tree-based group living. In 
addition, we did not find significant correlation between 
pair of individual genetic distance and geographic distance 
going beyond trees, suggesting group living on trees might 
not elevate inbreeding risk. 
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1. Introduction 


With the fragmentation of critical resources (eg. habitat or 
foraging site), a number of threatened individuals choose to live 
together to ensure their fitness by resource share, parental care 
or predator avoidance (Støen et al, 2005; Mckinnon et al., 2006). 
Nevertheless, group living likely dramatically increases the 
likelihood of mating among relatives and aggravate inbreeding 
risk (Charlesworth and Charlesworth, 1987; Vanthournout et 
al, 2016), which contributes a lot in genome-wide homozygosity 
and expression of deleterious recessive alleles (dominance 
hypothesis) and/or loss of heterozygous advantage (over- 
dominance hypothesis), thereby has several adverse impacts on 
population maintaining and individual survival (Charlesworth 
and Charlesworth, 1987; Madsen et al., 1999; Keller and Waller, 
2002; Charlesworth and Willis, 2009). 

To alleviate the adverse impacts, we need to manipulate 
group living and reduce the inbreeding level. The genetic 
relatedness between individuals should be lower than 
inbreeding risk level (e.g. 0.25) to maintain a sustainable 
population structure. Nevertheless, the association between 
group living and inbreeding risk is largely unknown. One 
recent study around the highly endangered red-cockaded 
woodpecker (Picoides borealis) found that inbreeding risk varies 
inversely with female natal dispersal distance, territory overlap 
among individuals was found to increase with the proportion 
of inbreed individuals (Schiegg et al., 2006). Uncovering the 
relationship between group living and inbreeding risk would 
contribute a lot to the conservation in endangered species 
(Cai et al., 2016; Brzeski et al., 2014; Hu et al., 2017), because an 


increasing number of species are endangered and forced to live 
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in small area with the aggravation of habitat fragmentation. 

The pit viper (Gloydius shedaoensis) provides an excellent 
model system for testing the relationship between group living 
and inbreeding risk. The snakes are endangered and occur 
exclusively on the Shedao Island (meaning “snake island” in 
Chinese), off the coast of Liaodong Peninsula in the Bohai Sea. 
The snake island is approximately 0.73 km’ and 53 km away 
from the nearest continent (Li, 1995). The population size of 
the snake is estimated at 20 281 (Li et al, 2007). The primary 
food source for the snakes is migrating birds in April-June and 
August-October every year (Li, 1995), although there are two 
small mammal species inhabiting on the island (bat Pipistrellus 
abramus and mouse Rattus norvegicus, Li, 1995). Small size of the 
island, high population density of snakes, and strong seasonality 
in food availability, force many snakes (2-9) to live on the same 
tree. Individual snakes consistently reused the branch of trees 
and prey on migrating birds (Shine et al., 2002), they would also 
mate on branches or nearby trees during group living (Li, 1995), 
which likely increase the inbreeding risk if they are relatives. 

The purpose of the present study was to examine whether 
group living on trees increase the inbreeding risk in G. 
shedaoensis. In particular, we asked whether snakes perched 
on the same tree were kin-related. The pairwise relatedness 
between snakes in and out of the same tree was estimated 
using microsatellite DNA phenotypes. We also compared the 
group construction based on different kin scenarios and trees, 
and detected the relationship between genetic distance and 
geographic distances to demonstrate the genetic relationships 
among individuals both in and out of trees. 
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2. Materials and Methods 


2.1 Field sampling To examine whether group living on trees 
increases the inbreeding risk in Shedao pit vipers (G. shedaoensis), 
a total of 45 snakes (15 males, 30 females from 10 trees) were 
captured using snake clamp from ten trees on Snake Islandin 
September, 2014 (Figure 1). After capturing, the snout-vent 
length (SVL, to 1 mm) and body mass (to 0.01 g) were measured 
using plastic ruler and electric scale, respectively, and the sex 
was determined by checking for the hemi-penis bulge. All 
snakes were permanently marked using PIT tags (Guangzhou 
HongTeng Barcode Technology Ltd., China), and a small 
drop of blood was collected from caudal vein using FTA card 
(Whatman International Ltd., UK) for subsequent genotyping. 
After carrying out all above procedures, snakes were released 
back to their dwelling tree. The coordinate of each capture site 
was recorded using GPS, and detailed sampling information is 
shown in Table 1. 


2.2 Microsatellite DNA collection We used microsatellite 
genotypes to estimate relatedness between pair of snakes. 
Because of lack in microsatellite markers specifically for G. 
shedaoensis, we first obtained eight microsatellite markers 
development for other Gloydius species by Goldberg et al. 
(2003), Holycross et al. (2002), Munguia-Vega et al. (2009) and 
Oyler-Mccance et al. (2005). The detailed information for those 
markers is listed in Table SL 

To increase the estimation accuracy, we specifically 
developed more microsatellite markers for G. shedaoensis. We 


Figure 1 a) an aerial view of Shedao Island; b) three Gloydius shedaoensis lived on the same tree; c) spatial positions of trees No. 1-9 where snake sam- 


pling was collected. 
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Table 1 Sampling information in Gloydius shedaoensis. 


Number of snakes 


Tree ID Latitude Longitude 
male female sexunknown sum 

No. 1 0 1 0 1 3857.006° N 12059.023^ E 
No. 2 4 4 0 8 3857.017°N 12059.033° E 
No. 3 3 9 0 12  3857.003' N 12059.008° E 
No.4 0 3 0 3 3856.999' N 12058.979° E 
No.5 D 2 0 4 . 3856.994 N 12058.989° E 
No. 6 1 3 0 4 3856.999° N 12058.980° E 
No. 7 p jl 0 3  3857.010 ' N 12059.021 E 
No.8 2 2 0 4 3857.039" N 12059.023° E 
No.9 1 p 0 3  3857.0388 N 12058.988° E 
No. 10 0 2 1 3 unknown 


acquired the transcriptome sequence of one G. shedaoensis from 
another project, and 40 microsatellite loci were identified and 
primers were designed to amplify those microsatellite DNA 
regions. Detailed isolation process has been described in Nie et 
al. (2015) and Li et al. (2015). Briefly, standard PCR amplification 
with optimized annealing temperature was conducted. PCR 
products were visualized on 8% denaturing polyacrylamide 
gels with silver staining, and the polymorphism of each locus 
was assessed. After multiple tests, nine loci were consistently 
amplifiable and polymorphic. Cross-species amplification was 
tested using three related Gloydius species, G. intermedius Strauch, 
G. strauchi and G. brevicaudus Stejneger. The detailed information 
for those newly developed markers is listed in Table 2. 

The genotypes of 45 individual snakes were collected using 
17 microsatellite DNA loci. In particular, we took a sample disc 
from the blood card for each snake and incubated it in NaOH 
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solution (20 mmol/L) for 25min under room temperature. 
The sample disc was then washed with TE-1 Buffer (10mM 
Tris-HCl, 0.1 mmol/L EDTA, pH 8.0) twice before being used 
in PCRs. PCRs were conducted with the same parameters as 
above. The forward primers were labeled with fluorescence 
dye of FAM, TAMRA, or HEX. PCR products were genotyped 
on an ABI 3730 genetic analyser (Sangon Biotech Co. Ltd., 
Shanghai, China). All eletropherograms were visually inspected 
and analyzed using GENEMARKER version 16 (SoftGenetics, 
State College, PA). All samples were genotyped at least twice to 
get consistent results. 


2.3. Genetic analysis Before genetic analysis, null alleles were 
checked for each locus with MICRO-CHECKER Version 2.2.3 
(Van Oosterhout et al., 2004). Loci of “S19”, “S11”, “S10” and “S22” 
were excluded due to low polymorphism. Genetic diversity for 
each locus was assessed using number of alleles (A), observed 
heterozygosity (H,), and expected heterozygosity (H,) with 
GENEPOP version 4.5.1 (Rousset, 2008). Deviation from Hardy- 
Weinberg equilibrium and linkage equilibrium was analyzed 
with GENEPOP Version 4.5.1. 

To examine whether snakes perched on the same tree were 
closely related from genetics, we calculated the mean relatedness 
between pair of snakes within trees using KINGROUP v2 
(Konovalov, 2010). Male and female were analyzed as overall or 
separately to ascertain whether genetic relatedness varied with 
sex. 

We then sorted individuals into subgroups by evaluating 
alternative partitions according to different likelihood scenario 
of kin relatedness. The overall likelihood for a given partition 


Table 2 Characterization of nine novel microsatellite markers for Gloydius shedaoensis and cross-amplification of those markers in three oth- 
er Gloydius species. Ta, optimal annealing temperature; N,, number of alleles; P, tests for Hardy-Weinberg equilibrium; H,, observed heterozy- 
gosity; H,, expected heterozygosity; F, forward; R, reverse; 4|, succeeded in amplification; x, failed in amplification; GIN, Gloydius intermedius 


Strauchi, GST, Gloydius strauchi, GBR, Gloydius brevicaudus Stejneger. 


Locus Primer sequence (5-3) Repeat motif Ta (C) N, Size (bp) P H, H, GIN GST GBR 
F: GGTTGAGAAACACTGCCTTAGC 

E OE E (AATO); 55 4 158170 004 077 060 d 4l 4i 
F: CAAAGCTACCCGCTTACCTG 

S19 (ACC), 55 1 227-236 000 100 051 4 4 4 
R: TCAGGATAAGTTTCGGGTCG 
F: AACCGGGCACAATATTTACG 

$24 (AGG), 55 2 280-286 000 097 051 4i 4l 4i 
R: CACGTGCTAAGTGTGGCAGT 
F.GTTCTGGCAGGAGGCACTAC 

SI (AATG), 55 3 181-185 — 100 033 032 ul 4 al 
R: AGGACTGTCACTTCCCATCG 
F: CAGTGACAGTGAAATCCAGCA 

Sit (AACT), 55 1 289 2 000 000 x d 4i 
R:CTGCTTCCTTCCATTTACCTT 
F: GGACCTCATCAACCAGCTTC 

S29 (AGG), 55 2 208-20 100 0.08 0.08 ul al 4 
R: CTTCTCCTGCAGGGTTCAAG 
F: AGCCACTTTCTCAACAACCTTC 

80 — RTTAATTGTOGATTGGGAGCTG (eL $ : d E Lu Me x y J 
F: GACAACCGCACGCACAAG 

S14 (AGCO); 55 6 282-20 049 036 037 x al 4 
R: ATGCAGGGATTCGTAAGGTT 
F: CTCCACCCTTTGCTCTTGTC 

$22 (ACO), 55 1 214 E 000 000 dj d di 


R: CTAAGTGGCTGAAACAGGGC 
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was calculated from the pairwise likelihood. To create partition, 
we first arbitrarily assigned one individual to the first subgroup 
while other individuals were added to the existing subgroup, 
one snake at a time. The placement of each individual is 
evaluated by placing it into available subgroups in turn as well 
as into a new subgroup, and taking the product of likelihood 
for corresponding pedigree relationships. The addition order 
was random and decided by the descending ratio search 
algorithm. We assigned another two different individuals into 
the first subgroup to examine the possible consequences of 
arbitrary reallocation, although no difference was found. All 
genetic group analysis was carried out using KINGROUP v2 
(Konovalov, 2010). 

To go beyond trees, an isolation-by-distance (IBD) pattern 
was examined to show the correlation between individual 
genetic relatedness and location distance. We used ‘mantel’ 
function in R package ‘vegan’ (Oksanen et al, 2019). The genetic 
distance was represented by F,,/(I-F,,) values. All geographic 
distances were calculated using R package ‘geosphere’ (Hijmans, 
2019). The analysis was only applied to site 1-9, because the 
coordinates of site 10 was missing. A Mantel test with 1000 
permutations was used to detect significant correlation between 
the genetic distances and geographic distances. 


2.4. Ethical approval All applicable international, national, 
and/or institutional guidelines for the care and use of animals 
were strictly followed. All animal sample collection protocols 


complied with the current laws of China. All animal procedures 


0.8 
a) 


oS m 
NO BR O^ 


Coefficient of relatedness (R) 
o 


-0.4 
0.6 Tree 
b) c) 
g 05 
203 
& o2 T 
2$ 01 0.0292 
$7 0.0246 
= 0 
So1 1 3 7 E 
8-02 
o 
'8-03 
E 0.4 
o 
o 


Tree 


| 
o 
N 
Hot 
wi 
— dm 
N 
oe] 


Coefficient of relatedness (R) 


Xiaoping WANG et al. 
Group Living and 
Inbreeding Risk 


245 


performed in this research were in accordance with the ethical 
standards of Experimental Animal Ethics Committee of the 
Chengdu Institute of Biology, Chinese Academy of Sciences 
(protocol number 2017005). 


3. Results 


3.1. Newly developed microsatellite loci Nine microsatellite 
markers were newly developed for G. shedaoensis (Table 2). The 
number of alleles in newly developed loci ranged from 1 to 6, 
the observed heterozygosity ranged from 0 to 1.00, while the 
expected heterozygosity ranged from 0 to 0.60. Five of them 
were significantly deviated from Hardy-Weinberg equilibrium, 
and six of them were crossly amplifiable in three other Gloydius 


species. 


3.2. Relatedness between snakes Mean relatedness between 
pair of snakes for all individuals, males and females perched 
on the same tree were demonstrated in Figure 2. In high 
proportion of trees (8/10), snakes did not show close genetic 
relationships, except the tree No. 7 (r = 0.30 + 0.16, n = 3) and 
tree No. 10 (r = 0.35 + 0.26, n = 3). The mean relatedness between 
males or females perched on the same tree was similar with the 
overall analysis, except the tree No. 3 (r = 0.40 + 0.02, n = 3) for 
males. 


3.3. Kin group reconstruction Kin group reconstructions 
based on different relatedness likelihood were demonstrated 


in Figure 3. Under parental-offspring scenario, all snakes were 


mum the coefficient of relatedness 
between individuals perched on a 
tree 


—— the average coefficient of 
relatedness among the trees 


0.06439 
F 10 0.00986 . 
the average coefficient of 
relatedness among the whole 
population 
0.6 
0.4 
02 | -0.0053 
0 | B —0.0149 
8 
-02 
-0.4 
—0.6 
Tree 
-0.8 


Figure 2 Average coefficient of relatedness (R) between pair of a) all individual Gloydius shedaoensis perched on the same tree; b) males perched on the 
same tree; c) females perched on the same tree. The solid line represents the average coefficient of relatedness between pair of individuals among the 
trees, while the dotted line represents the average coefficient of relatedness between pair of individuals from all trees. 
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Snake 

Bi Tree No.2 
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Figure 3 Kin group reconstruction based on scenarios of parental-offspring, full-sibling, half-sibling and cousin. Snake groups based on trees were pro- 
vided as reference. The same color bar represents Gloydius shedaoensis from the same tree. Snake from tree No. 1 was not included because there is only 


one pit viper from this tree. 


sorted into 15 subgroups. To be conservative, we focused on the 
first five subgroups. Two snakes from tree No. 2 were sorted 
into first two subgroups respectively. Two snakes from tree 
No. 3 and one snake from tree No. 4 were sorted into the third 
subgroups. One snake from tree No. 3 was sorted into the fourth 
subgroup, while the fifth subgroup consisted of one snake from 
tree No. 3 and one snake from tree No. 6. 

Under full-sibling scenario, eleven subgroups were created. 
One snake from tree No. 4 was sorted into the first subgroup, 
one snake from tree No. 5 was sorted into the second subgroup. 
The third subgroup consisted of one snake from tree No. 3, 
one snake from tree No. 5 and one snake from tree No. 6. The 
fourth subgroup consisted of one snake from tree No. 2, one 
snake from tree No. 3 and one snake from tree No. 7, while the 
fifth subgroup consisted of one snake from tree No. 7. 

Under half-sibling scenario, nine subgroups were created. 
One snake from tree No. 2 was sorted into the first subgroup. 
The second subgroup consisted of one snake from tree No. 3, 
two snakes from tree No. 5 and one snake from tree No. 6. 
The third subgroup composed of one snake from tree No. 2, 
one snake from tree No. 3 and one snake from tree No. 7. The 
fourth subgroup consisted of one snake from tree No. 2, one 
snake from tree No. 3, two snakes from tree No. 5 and two 
snakes from tree No. 7, while the fifth subgroup was made up 
of one snake from tree No. 8. 

Under cousin scenario, six subgroups were created. One 
snake from tree No. 5 and one snake from tree No. 6 were 
sorted into the first subgroup. The second subgroup composed 
of two snakes from tree No. 2, one snake from tree No. 3, two 
snakes from tree No. 6 and two snakes from tree No. 7. One 
snake from tree No. 8 was sorted into the third subgroup. The 


fourth subgroup consisted of two snakes from tree No. 2, five 
snakes from tree No. 3, one snake from tree No. 5, one snake 
from tree No. 7 and two snakes from tree No. 8, while the fifth 
subgroup composed of one snake from tree No. 2, one snake 
from tree No. 3, two snakes from tree No. 9 and one snake from 
tree No. 10. 

No matter which scenario, the group construction 
was different with the tree-based group living, suggesting 


individuals living on the same tree was not kin-related. 


3.4. Isolation-by-distance The correlation between individual 
genetic distance and geographic distance was shown in 
Figure 4. We found no significant correlation between pair of 
individual genetic distance and geographic distance (r = 0.12, P = 
0.34). 


0.15 P=0.34; r=0.12 


0.10 


FJ 0-Fy) 


0.05 


0.00 


0 30 60 90 
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Figure 4 Isolation-by-distance pattern between pair of individu- 


al Gloydius shedaoensis genetic distance and geographic distance. 


| No. 2 


4. Discussion 


Using genotypes from 13 microsatellites DNA markers (eight 
from other Gloydius species, and five from newly developed 
microsatellites for G. shedaoensis), we found snakes living on the 
same tree were largely not kin-related. We have three evidences 
for this. First, except the tree No. 7 and No. 10, snakes on the 
same tree did not show close pairwise genetic relationships. 
Second, the group constructions based on parental-offspring 
scenario, full-sibling scenario, half-sibling scenario as well as 
cousin scenario, were not consistent with snakes group based 
on trees. Third, genetic distance and geographic distance was 
not correlated between pair of individuals. This is very different 
to social group living lizards of genus Egernia, among which 
23 of 30 reported species are found to aggregate based on kin 
relatedness (Chapple, 2003), and E. striolata was found to form 
a tree-dwelling population in group of close kin (Duckett et al., 
2012). 

None kin-related group living in G. shedaoensis implies that 
snakes might have taken crucial ways in avoiding inbreeding 
risk before tree occupation. Offspring and adult dispersal after 
hibernation likely play important roles in inbreeding avoidance. 
As an island-dwelling species, snakes G. shedaoensis confront 
many challenges from small living area and high population 
density. Adult G. shedaoensis choose to stay in most of time to 
avoid social conflict. The average daily movement distance of 
male and female is 1.76 m and 1.99 m, respectively, and a radio- 
telemetric study shows that adult snakes generally remain 
within 50 m of their original site, with no tendency to move 
further (Shine et al, 2003). Nevertheless, offspring G. shedaoensis 
was found to disperse intensively during breeding season, out 
of 11 marked offspring, six of them have low recapture rate. In 
addition, adult G. shedaoensis are very active after hibernation, 13 
of followed adults have average of 4.59 + 5.48 m (range from 0 
to 39 m) in dispersal distance per day (Yayong Wu, unpublished 
data). Dispersal behavior has been found to contribute a lot 
in inbreeding avoidance in many other species (Perrins and 
Goudet, 2001; Olsson and Shine, 2003; Le Galliard et al., 2006; 
Schiegg et al, 2006). For example, in endangered red-cockaded 
woodpacker P. borealis, natal dispersal distance is negatively 
associated with inbreeding risk, suggesting natal dispersal could 
alleviate inbreeding pressure (Schiegg et al., 2006). 

Alternatively, female G. shedaoensis likely avoids inbreeding 
through strict mate choice. This is inferred from the interesting 
courtship behavior described by Li (1995). Before copulation, 
male G. shedaoensis always tend to conquer females using 
forcible ways, such as twining and biting, while females resist 
males using the same way. As soon as copulation starting, 
females frequently nod their head (Li, 1995). We assumed that 
head nodding behavior in female G. shedaoensis likely functions 
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as important signal of female acceptance because this behavior 
only happened during copulation. We still need more direct 
experiment to verify this, because female mate choice is very 
tricky in reptiles, although there is one evidence on cryptic mate 
choice in Swedish sand lizards (Lacerta agilis) (Olsson et al, 1996). 
We admitted that the relationship between group living and 
inbreed risk we obtained in G. shedaoensis should be cautious, 
because in most trees, the number of snakes sampled was low, 
which likely has consequences on the accuracy of relationship 
estimation. We need to add more snakes from more trees to 
verify this in the future. 


Conclusion 


Our study evidenced that tree-based group living in adults G. 
shedaoensis did not aggravate inbreeding risk, because snakes on 
the same tree are not kin related. This implies that snakes may 
have taken crucial ways in avoiding mating with relatives, such 
as offspring natal dispersal, adult dispersal after hibernation 
and strict mate choice, although much of those are speculated 


and warrants careful future research. 
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Appendix 


Table S1 Information of eight microsatellite loci from published literature. 
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Locus Repeat motif Primer sequence (5’-3’) Ta(C) Allelic Size (bp) Fluorescence labeling Reference resource 

Crti37.— (GT), (GA), ue c HUMUM mo 230-246 FAM Munguia-Vega et al. (2009) 
MFRD5 (TG)23 x Prosa etn ied 60 168-174 HEX Oyler-Mccance et al. (2005) 
Crti95 (CA)22 es = 2 pce. ae ee 50 164-180 FAM Munguia-Vega et al. (2009) 
Ch7-87 (CA)12 Beyer E dies perenne E G 60 146-162 HEX Villarreal et al. (1996) 
Crtil2A — (CA22 es i TUR em = * mi. G 50 210-214 TEMRA Munguia-Vega et al. (2009) 
Ch5A (CA)17 : d uU A 60 152-166 HEX Villarreal et al. (1996) 
Crtil2 (CA)14 x E CORE RUN 54 228-242 TEMRA Goldberg et al. (2003) 
CwB6 (GA)19 FOTCITTTACGCCCACCACTTTA; 56 106-110 FAM Holycross et al. (2002) 


O 


R: CCCCGCTAACCTTTGCTCAG 
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